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We present three-body nonrelativistic calculations within the framework of a potential model for 
the kaonic cluster K~pp using two completely different methods: the method of hyperspherical 
harmonics in the momentum representation and the method of Faddeev equations in conhguration 
space. To perform a numerical benchmark, different NN and antikaon-nucleon interactions are 
applied. The results of the calculations for the ground state energy for the K~pp system obtained 
by both methods are in reasonable agreement. Although the ground state energy is not sensitive 
to the NN interaction, it shows very strong dependence on the KN potential. We show that the 
dominant clustering of the K~pp system in the configuration A(1405) -I- p allows us to calculate 
the binding energy to good accuracy within a simple cluster approach for the differential Faddeev 
equations. The theoretical discrepancies in the binding energy and width for the K~pp system 
related to the different NN and KN interactions are addressed. 


I. INTRODUCTION 


The theoretical and experimental studies of composite systems of mesons and baryons is a challenging issue in 
hadronic and nuclear physics. The study of the KN interaction and the formation of the A(140M as a quasi-bound 
state of KN started in the 60s. Early suggestions for kaon condensate in dense matter motivated the 

search for bound states of kaons in nuclei, since the kaon-nucleus interaction could answer the question of whether 
kaon condensation takes place in the core of neutron stars. Kaonic systems have attracted much attention in the 
last decade. In Ref. Q it was shown that kaonic nuclei could be produced by the {K~, p) and {K~, n) reactions. 
However, the renaissance in this field started after the prediction of deeply bound K— states by calculations performed 
by Akaishi and Yamazaki in Ref. Q and followed by several subsequent publications imi. The experimental search 
was initiated at the KEK 12 GeV proton synchrotron in Japan by the finding of peaks in the nucleon spectra of 
A'“-absorption in "^He 0[ll. 

The KN interaction has been found to be strongly attractive in the isospin 7 = 0 channel and, as a consequence of 
that, nucleons are attracted by the K meson to form dense nuclear states. In 1963, in Ref. [l^ Nogami first examined 
three possible isospin configurations of the KNN system and discussed the possible existence of the bound state using 
a rather crude calculation. The primary focus was on the lightest kaonic system K~pp because it has two K~p pairs 
with isospin 1 = 0 that can be considered as a building block of a three-body kaonic cluster. Regardless of the fact 
that two protons are not bound, there is an expectation that KN interactions can be strong enough to bound the 
K~pp system. Interestingly enough, even two protons can be bound by a single K meson and it is assumed that this 
would be the strongest binding cluster among the three-body systems! 

On the experimental side, several experiments have been performed to search for the kaonic K~pp cluster using 
various reactions. The first measurement reported by the FINUDA collaboration was the measurement of the invariant 
mass distribution of the p produced by AT “—absorption by ®Li, ^Li, and at the e“'"e“ rings of the DA^NE machine 
in Frascati da. The analysis of the invariant mass distribution of the Kp produced by AT — absorption gave the 
value (2255±9) MeV/c^ for the mass of the peak, corresponding to the binding energy 115t5(stat)t4(syst) MeV 
and the width 67l4^(stat)l3(syst) MeV. However, the authors of Refs, [l^, [l^ claimed that the FINUDA data 
can be explained without postulating the existence of the K~pp bound state and suggested that the peaks of the 
proton spectra come from the K~ absorption on a pair of nucleons, leaving the rest of the nucleons as spectators. 
Interestingl y en ough, even though such a mechanism cannot be completely excluded, there exists the further counter¬ 
arguments against the interpretations of stopped-AT” experimental data of the FINUDA group by [l^. The 
OBELIX experiment at LEAR-CERN also suggested the observation of the K~pp state with invariant-mass 
spectroscopy of stopped p reactions on ''^He. However, a less significant signal was reported in the stopped-p reaction. 
After the FINUDA collaboration observation [13: there have been reports of possible experimental evidence for the 
kaonic K~pp state in heavy ion collisions. Results were obtained by the FOPI group experiment [2^ at GSI in the 
study of Ni-|-Ni and Al-I-Al collisions with the production of Np and Ad. The DISTO collaboration at the SATURNE 
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machine have reanalyzed their dataset of the experiment on the exclusive pp —>■ pAK~^ reaction at 2.85 GeV to search 
for a strongly-bound compact K~pp state formed from the p-|-p —>■ K~^ + {K~pp) reaction [2^, The enormously 
large bump structure observed in the K'^ missing mass and the pA invariant-mass spectra at 2.85 GeV indicate a 
possible candidate for the formation of a compact K~pp cluster with large binding energy. The binding energy and 
width were determined to be 103±3(stat)±5(syst) MeV and 118±8(stat)±10(syst) MeV, respectively. Therefore, two 
collaborations da and [ 2 ^, [ 2 ^ claim that K pp is strongly bound with a binding energy of more than 100 MeV, 
although the decay width is rather different between these two. The first search results using the yd —^ A"+7r“X 
reaction at the range of photon energy 1.5-2.4 GeV were reported in Ref. [ 2 ^. A statistically significant peak structure 
was not observed in the differential cross section of the AT+tt” photo-production in the region from 2.22 to 2.36 GeV/c^ 
in the inclusive missing mass spectrum. However, the upper limits of the differential cross section of the K~pp bound 
state production were determined with the assumed widths of 20 MeV, 60 MeV and 100 MeV. 

Ref. [ 2 ^ reported on preliminary results of the J-PARC E15 experiment aiming to search for the bound state of 
K~pp via the in-flight ^He(A'”, n) reaction at 1 GeV/c and the first data collection was performed with 5x10® incident 
kaons on the ®He target. The experiment investigates the K~pp bound state exclusively both in the formation via 
missing-mass spectroscopy and its decay via invariant-mass spectroscopy using the emitted neutron and the expected 
decay, K~pp —> A pp —> 'K~pp, res pec tively. Data analyses of the semi-inclusive and exclusive channels are in 
progress. At the same time, in Ref. [27j| have been reported results for the in-flight kaon-induced reaction on ^He: 
K~ -|-^He—^ K~pp + n. In particular, the semi-inclusive analysis of the neutron missing-mass spectrum, no significant 
peak was observed in the region corresponding to K~pp binding energy larger than 80 MeV, where a bump structure 
was reported in the Ap final state in different reactions by FINUDA and DISTO collaborations. Thus, the situation 
is still controversial and the existence of the K~pp bound state has not been established yet. New experiments using 
different reactions could help to resolve this controversial situation. 

Such experiments are being planned and performed by the HADES collaboration [^ . It was recently proposed 
to search for kaonic nuclei at the future SuperB factory, which will be built in the Tor Vergata University Gampus 
in Rome. The K~pp cluster can be identified through its Ap decay mode, searching for a narrow peak at a mass of 
about 2.25 GeV/c® in the Ap invariant mass spectrum and measuring the A — p angular correlations, which can give 
important hints on the nature of the event [^. One of the advantages of searching for light A'”—nuclear clusters at 
the future SuperB factory is that the search could be extended from the nuclear medium to the vacuum, looking for 
its production in the strong decays of T(1S) by taking advantage of the high luminosity of these machines. Also, all 
problems related to the influence of the medium as well as the hnal state interaction will hopefully be avoided. 

On the theoretical side, there are many reports related to K —nuclear systems. Today, refinements of the existing 
methods for studying few-body systems, on one hand, and developments of new methods, on the other hand, and 
advances in computational facilities enable very precise calculations for few-body A'—nuclear systems. The light 
kaonic KNN cluster represents a three-bo^ system and has been treated in the framework of various theoretical 
approaches such as variational methods (ol-fisl. Isil - ls^ , the method of Faddeev equations (sol - ls^ and the method of 
hyperspherical harmonics (HH) ■ 

In the framework of the Brueckner-Hartree-Fock theory, by defining the KN p-matrix in a nuclear medium as well 
as using antisymmetrized molecular dynamics for the spatial functions, the K~pp cluster has been studied in Refs. 
i-H, Variational calculations of the AT pp system using several versions of energy-dependent effective KN 

interactions derived from chiral SU(3) dynamics were performed in Refs. [s^. A different ansatz was used in Ref. 

for the trial wave function of strongly correlated nucleons in the variational calculations of the binding energy and 
width for the K~pp cluster. Recently in Ref. the K~pp cluster was investigated by combining a coupled-channel 
Gomplex Scaling Method with the Feshbach resonance method allowing one to effectively treat a coupled-channel 
problem as a single channel problem. The authors use an energy-dependent chiral-theory based potential and employ 
the correlated Gaussian function [38l |. They report that the K~pp cluster is shallowly bound with a binding energy 
of around 20-35 MeV and the half decay width ranges from 20 to 65 MeV. 

The strongly bound nature of the K~pp system is shown in Refs. by using the coupled-channel Faddeev 

equations in the Alt-Grassberger-Sandhas form and solving for all two-body separable potentials to find the 
resonance energies and widths. Similar calculations have been done in Refs, [dll, and [d^. The K~pp cluster is 
studied in Refs. [dB - ldSj l using the fixed center approximation to the Faddeev equations that relies upon on shell two 
body amplitudes. Thus, the structure and the production mechanism of the K~pp bound state have been investigated 
using various theoretical approaches. All the aforementioned approaches predict the existence of a bound state for 
the KNN system. While on the one hand, various theoretical approaches predict a strong bound state of the K~pp 
system, on the other hand, various theories based on chiral dynamics predict a shallow bound state for the K~pp 
cluster. The predicted values for the binding energy and the width are in considerable disagreement: 9-95 MeV and 
20-110 MeV, respectively. The serious theoretical discrepancies concerning K~pp binding and width essentially come 
from different starting ansatzes: phenomenological AT-nucleon potential constructed using KN scattering data and 
the data of kaonic hydrogen atomic shift together with the ansatz that A(1405) resonance is a bound state of KN 
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system, a chiral ansatz for KN interaction that leads to a substantially shallow KN potential, and the calculation 
methods. Further theoretical investigations are apparently needed to resolve this controversial situation that perhaps 
stems from the ambiguity of the KN interaction, the importance oi NN interaction or the different procedures for 
three-body calculations. 

We present three-body nonrelativistic calculations within the framework of a potential model for the kaonic cluster 
K~pp applying two completely different methods: the method of hyperspherical harmonics in momentum represen¬ 
tation and the Faddeev equations in configuration space. Calculations for a binding energy and width of the kaonic 
three-body system are performed using three different potentials for the NN interaction, as well as two different 
potentials for the description of the kaon-nucleon interaction. One is the energy independent phenomenological KN 
potential from Ref. and the other is the energy dependent chiral KN interaction. Such approach allows one to un¬ 
derstand the dependence of the bound state and width of the kaonic three-body system on the method of calculations, 
the importance of nucleon-nucleon interaction, and the key role of the kaon-nucleon interactions. 

This paper is organized as follows. In Sec. II we describe the effective potentials of the K~p interaction and 
present the theoretical framework to solve the three-body problem for the K~pp system. We apply the method of 
hyperspherical functions in the momentum representation for solving the three-body Schrodinger equation and the 
method of Faddeev equations in configuration space to find bound states for the kaonic K~pp system. The results of 
the three-body calculation for the binding energy and the width for the K~pp cluster and the analysis of the structure 
of this system are presented in Sec. 3. We also compare our results with those obtained within variational methods, 
the method of hyperspherical functions in the coordinate representation and the Faddeev equations in momentum 
representation. Finally, in the Sec. 4 we present the summary of this work and draw conclusions. 


II. THEORETICAL FRAMEWORK 
A. Interactions 

The Hamiltonian of the three nonrelativistic particles for the KNN system reads 

H = T + VnN + T/caTi + (1) 

where T is the operator of the kinetic energy, Vnn is the nucleon-nucleon potential and + ^kn 2 

pairwise effective antikaon interaction with the first and second nucleon, respectively. The effective interactions of the 
KN, KN, KK and KK two-body subsystems are discussed in detail in Refs. (sl. 1^. IT^. 1^. 1^. [134^ . Below, we use 
two effective KN interactions that were derived in different ways. The effective KN interactions can be constructed 
based on a phenomenological approach so as to reproduce the existing experimental data for the KN scattering 
length, the mass and width of the A(I405) hyperon and the Is level shift caused by the strong KN interaction in the 
kaonic hydrogen atom. The effective KN interaction can be derived within the chiral SU(3) effective field theory, that 
presents the low-energy realization of QCD with strange quarks and identifies the Tomozawa-Weinberg terms as the 
main contribution to the low-energy KN interaction [^. The potential for the description of the KN interactions 
was derived in Refs. iEl phenomenologically using KN scattering and kaonic hydrogen data and reproducing the 
A(1405) resonance as a K~p bound state at 1405 MeV. We refer to this as the Akaishi-Yamazaki (AY) potential. 
The AY potential is energy independent. The other KN interaction given in Ref. [s^ was derived based on the 
chiral unitary approach for the s—wave scattering amplitude with strangeness S' = — 1 and reproduces the total cross 
sections for the elastic and inelastic K~p scattering, threshold branching ratios, and the ttE mass spectrum associated 
with the A(1405). We refer to this energy dependent potential for the parametrization as the HW potential and 
in calculations with this potential, we are following a procedure [§3 . . Both potentials are local and constructed in 

coordinate space. 

The potentials for both above mentioned effective KN interactions can be written in the one-range Gaussian form 
as 


yRNir) = ^ U^exp 

7 = 0,1 

where r is the distance between the antikaon and the nucleon, b is the range parameter and Pj^j^ is the isospin 
projection operator. The values of the potential depth and for each interaction are given in Refs. [8l.fl3ll33| 
and the range parameter is chosen to be 5 = 0.66 fm for the AY potential and b = 0.47 fm for the HW potential. 


- {r/by 


KN^ 


( 2 ) 
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To describe a nucleon-nucleon interaction, we use three NN potentials: the realistic Argonne V14 (AV14) [^, the 
semi-realistic Malfliet and Tjon MT-I-III (MT) potential and the Tamagaki GSRS potential [blj that we refer to 
as the T potential. The T potential was used in Ref. m for variational calculations and that allows us to compare 
our results with the latter work. We chose the MT s-wave potential to compare how the binding energy depends 
on the shape of the s-wave NN interaction. The MT s-wave potential has only two components in spin-isospin 
channels (s,t) = (l,0), and (s,t)= (0,1). The potential allows one to reproduce and obtain an acceptable description 
of experimental data for bound states and scattering in three and four nucleon systems The realistic AV14 

potential perfectly fits both the pp data, as well as the np data, low-energy nn scattering parameters and the deuteron 
properties. This potential was chosen in order to study the influence of the orbital partial wave components of the 
NN interaction (by the non central operator) with I > 0 on the energy of the K~pp cluster. In our considerations, 
the other components (spin-orbital, tensor and other) of the realistic potential have not been taken into account 
due to their weak effect on the system. Moreover, in the model with one spin-isospin channel (s,t)=(0,l), the other 
components of the potential do not make contributions to the nucleon-nucleon interaction. Thus, the use of all of 
these potentials allows the validity test against various NN potentials. 


B. Formalism of HH in momentum representation 

Let’s introduce the trees of Jacobi coordinates in configuration and momentum spaces for a system of three particles 
with unequal masses mi, m 2 , and m 3 having positions ri, r 2 and ra and momenta ki, k 2 , and ks as follows 


mjmk 
rrij + mk 


(rj - Tfc), 


yz = 


l mj{mj + mk) ( mjrj -|- m^rfe \ 

V M V * mj + mk ) ' 


R = (miri -I- m 2 r 2 -l- msrs), M = mi -|- m 2 = m 3 , i ^ j ^ fc = 1, 2, 3. 
The conjugate sets of the Jacobi momenta (qi,Pi,P) for the partition i are defined as: 


(3) 


Qi = 


1 


= (mfekj - mjkfc). 


^mjmk{mj + mk) 

Pi =- ! ^=[-m,(kj- -I- kfc) -I- (mi -I- mfc)ki]. 


^Jm^M(mi -I- mk 


P, = 


1 




(ki -t- k 2 -k ka), i ^ j ^ fc = 1, 2, 3. 


(4) 


We introduce the hyperspherical coordinates as the hyperradii g and >«■, and two sets of hve angles denoted by Oi = 

(ai,Xi,yi) and fli = (/3i,qi,Pi) which define the direction of the vector p and vector , in the six dimensional 
configuration and momentum spaces, correspondingly, so that 


P = + yh = /ocosoi, Pi = psinai, 

dxidyi = p^dpsin^ acos^ adad%dyi = p^dpdVLp, 
^ = \/pi + dh Pi = >^sin/3i, Qi = xcos/di, 
dpidqi = p^dp sir? (3 cos^ I5dl3dpidx\i = ^c^d^cdfli. 


(5) 

( 6 ) 


The hyperradii g and >c are invariant under three-dimensional rotations and independent of the partition i. 

One can write the Schrodinger integral equation describing three bound particles in the momentum representation 
as 


^(p,q) = - 


1 


(27r)e 


G(p,q)<pq |4h23|pq>4'(p ,q)dp'dq'. 


(7) 
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where 


< p q 1^4231 pq > 


1 


Vi 23 exp 


*(q- q )x+*(p 



dxdy, 


( 8 ) 


is the Fourier transformation of the V 123 , which is defined as the sum of the pair-wise nucleon-nucleon and effective 
antikaon-nucleon interactions. In Eq. 0 , the Green function has the form 


G(p,q) 


2M 1 
It? + Kq' 


(9) 


where and E is the bounding energy of the three-particle system. 

For the next step, we follow the procedure given in Refs. [651 166| and expand the wave function of three bound 
particles in terms of the antisymmetrized hyperspherical harmonics cr, r) in the momentum representation: 

o, r), ( 10 ) 

f-llp Iq 


where /r is the grand angular momentum, L is total orbital momentum and Ip and Iq are the angular mo¬ 
menta corresponding to the momenta p and q. In Eq. (nni) the antisymmetrized hyperspherical functions 
cr, r) are written as a sum of products of spin and isospin functions and hyperspherical functions 
(^Ipl^rnpUiql LM) where are hyperspherical harmonics, using the 

rripTriq 

Raynal-Revai coefficients [g^. The HH are the eigenfunctions of the angular part of the six-dimensional 

Laplace operator in momentum space with eigenvalue fj,{p + 4), and are expressible in terms of spherical harmonics 
and Jacobi polynomials [6^[63|. By substituting Eq. (fTTO into the integral Schrbdinger equation in the momentum 
representation © and taking into account © one obtains a system of coupled integral equations for the hyperradial 
functions : 


^ (^plql^plq) j pdpJp+2{^p)Jp'P)^ 

(aj VA{xi)Yf 


( 11 ) 


where the index A is related to the type of interaction A G NN, KN. In Eq. yjj ^Q-p^ql^p^qYpL are the Raynal-Revai 
coefficients [g^ for the unitary transformation of HH from one set of Jacobi coordinate to another, J^+ 2 (>!rp) are the 
spherical Bessel functions and 


= fV'p'''(cosai)''’(sinaj)'5P,(p+^/^l«+^/^(cos2Q;i), 


_ / 2n!(^-I-2)(n-I-/p-I-/g-I-1)! 

~ YT{n + h + i/2)T{n + l2 + il2Y 


( 12 ) 


where 2n = p — Iq — Ip, n=l, 2, 3,... and Pn^'' is the Jacobi polynomial. By solving the coupled integral equations 
(CD one can find the hyperradial functions for a given L. 


C. Formalism of Faddeev equation in configuration space 

The wave function of the three-body system can be obtained by solving the Schrbdinger equation with the single 
channel Hamiltonian ([T} . Alternatively, in the Faddeev method the total wave function can be decomposed into three 
components: = $1 J- 4>2 J- $3 [^. The Faddeev components 4*^ correspond to the separation of particles into 

configurations i J- (A:/), i fc 1 = 1, 2,3. Each Eaddeev component = <I>i(xi,yi) depends on its own set of the 
Jacobi coordinates The components satisfy the Eaddeev equations in the coordinate representation written in 
the form [toI - It^ : 

{Hq + v^Xi) - E) $*(xi, yi) = -vYm) (4’fc(xk, yk) + 4’z(xi, yi)), (13) 
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where i k ^ I = 1,2,3. Hq = — ^(Ax^ + AyJ is the kinetic energy operator, Mi is the reduced mass for 
the corresponding Jacobi coordinates dS]) and vi is the potential acting between the particles {kl), i ^ k ^ 1. The 
orthogonal transformation between three different sets of the Jacobi coordinates has the form: 


X* 




Cik Sik 

^ik ^ik 



^ik + ^ik — 1 : 


(14) 


where 


Cik = — 


rrnmk 


{M - mi){M - nik) 


S^k = (- 1 )'^ *sign(fc - i)^l- 


The bound state of the K pp system can be obtained by solving the Faddeev equations in configuration space m- 
In this case the total wave function of the K pp system is decomposed into the sum of the Faddeev components 
U and W corresponding to the {pp)K~ and {K~p)p types of rearrangements: 41 = U + W — PW, where P is the 
permutation operator for two identical fermions. For a three-body system that includes two identical fermions the set 
of the Faddeev equations (IT^ can be reduced to the system of two equations for the components U and W [zl, [t^ : 


(iJo“ + Vpp - E)U = -Vpp{W - PW), 

{H^ + Vk-p - E)W = -Vk-p{U - PW), 


where the potentials for pp and K p pairs are defined by Vpp and Vk-p, respectively. The partial wave analysis of 
the differential Faddeev equations (DFE) (fT51) in the LS basis is performed by the general scheme described in Refs. 
[ 2 l[zl,[ 2 l. The LS basis allows us to restrict the model space to the states with total spin S = 0 (when the spin 
projections of protons are anti-parallel) and total isospin I = \. The possible isospin configurations with J = 1 are 
not taken into account in our calculations. According to the evaluations of different authors the total contribution 
of the 1=1 configuration is about 5% [l^. We consider two cases for the solution of the system (flSl) : the s-wave 
approach and the cluster approach. In the latter case we assume the dominant clustering of the K~pp system in 
the form A(1405) +p, that allows us to present the Faddeev components as the product of the eigenfunctions of the 
two-particle subsystems and functions of the relative motion of pair and the third particle. 


1. s-wave approach 

Let us consider the s-wave approach. Above, we mentioned that the nucleon-nucleon MT and T potentials and 
K~p interactions applied have only the s-wave components for the spin-isospin state s=0, t=l in the pair NN and 
t=0,l in the pair K~p. We consider the s-wave approach for the differential Faddeev equations (fTSl) . The s-wave 
DFE can be written as 




Vppix) - E)U{x,y) = -vl 


(x) I du^^AW{x[,y[), 


-1 




{HS^^VK-p{x)-E)W{x,y) 



■ 1 

/ du-^A^U{x'2,y'2)+ 
J x'2y2 
-1 


(16) 


1 

+ / du^^DW{x2,y2) , 

J X 2 V 2 

-1 

where the numbers of particles in the system are labeled as 2 and 3 for protons and 1 for the kaon, m 2 = m 3 = m 
is the nucleon mass, and mi is the mass of the kaon. Vpp(x) is the singlet-triplet (s =0 and t=l) component of the 
NN potential, Vk-p =diag{i;J-_p, with and the singlet and triplet isospin components of KN 










7 


potential, respectively. In Eqs. (IT6l) u=cos(xy), where xy is the angle between x and y, and we define 




TTU _ _ 

° 2 m(i) 


-d?.- 






-dl- 






D = 


TTW _ _ 

° “ 2M(2)^^ 2m(2) 

1 ^ \ 

2 2 

^ _i r 

2 2 / 




W{x,y) = {W{x,y),W\x,y))'^. 

In (fT6)l . the appropriate transformation of coordinates and reduced masses are given by 


^ ((^ 
\ m 


\2 2 
c)2 + y2_p2 


m 


\ 1/2 

\ , rrii // ™ 2 r. Ill \ 


1/2 


(17) 


(18) 


.m 


( 2 ) 


* 2 = {^xr+y^ + 2- 


7,(2) 


m 


-xyu 


1/2 


2 / 2=2 


M(2) 


r)2+y2-2 


m 




\ 1/2 


= ( 


where to(^) = 


7.(2) 


mi 


c)" + y "-2 


m 


( 2 ) 


mi 


-xyu 


1/2 


1/2 = 


mi 


('■ 


mi 


M(b= 2 ^^, 

mi H- 2m 


M( 2 ) = 


M( 2 ) 
m(mi 


x )2 + y" + 2 


mi n1/2 


+ m) 


^( 2 ) _ mim 
mi+m 


mi + 2m 

are reduced masses the pp pair and the pair and third particle {K~ 


Spin-isospin 


configurations of the subsystems of the K pp system, corresponding to (fT6)l are graphically presented in Fig. [T] 


K O 


P 

a) 




FIG. 1: Spin-isospin configurations of the subsystems of the K pp system, a) K -\-{pp) with the singlet spin state of the pair 
{pp). b) {K~p) -I-p with the singlet isospin state of the pair {K~p). c) {K~p) p with the triplet isospin state of the pair 
(K-p). 


2. Averaged potential 

The K meson combines two nucleons into the bound state for two isospin configurations that are energetically 
favorable. The effective KN interactions have strong attraction for the 1 = 0 channel and weak attraction in the 
7 = 1 channel. In this section, we consider an approach for the case when one of the pair potentials has different 
components for spin-isospin states. Following, for instance Refs. [lH and [b^, for Fqs. (HU) and dm), one can consider 
an effective potential obtained by averaging the initial potential over isospin variables. Fach component of the initial 
potential is substituted by the averaged potential. The isospin averaged potential Vx-pir) is defined by the following 
form: 

VK~p{r) = \v^£°p{r) + ^v^£}p{r). (19) 

This potential has moderate attraction in comparison with the component for 1 = 0. Using this potential we reduce, 
for example, the set of equations (HH) to two equations. The averaged potential replaces the Vx-p so that 

Vx-p = 

and W{x,y) becomes W{x,y) = —^lUi — ^W 2 . Note that this simplification changes the two-body threshold in Eq. 
(US), which is not related to the Vx-p bound state as A(1405). The value of the two-body threshold for the effective 
potential is about -9.6 MeV instead of -30.2 MeV for the AY potential (7 = 0). 























3. Cluster approach 


Within the formalism of the Faddeev equations, one may use dominant clustering of the K~pp system in the form 
A(1405) + p to calculate the binding energy. In the framework of such an approach, the s-wave Faddeev components 
are decomposed as a product of the eigenfunctions of the Hamiltonians of two-particle subsystems and functions of 
the relative motion of the pair and the third particle. 

Here, we present a brief description of the method based on the general form for the Faddeev equations (fT^ . The 
Faddeev components Lion a=l,2,3 are written in the following form [t^ : 

Ua{x,y) = 4'a{x)fa{y). ( 20 ) 

In (I20L the functions (j)a are the solutions of the two-body Schrodinger equations for subsystems with minimal 
eigenvalue 


n? \ 

-^-^dl+Va{x)\ (j)a{x) = £ra(/)a(x). 

Substituting (1201) into the Faddeev equations and projecting, one may obtain the set of integro-differential equations 
for the functions fa{y) describing the relative motion of clusters: 




2Mh 


fa{y) 


= Va{x) J du|-^((>/3(x^)//3(2/^) + | ^ 


( 21 ) 


In these equations, (.|.) means the integration over the variables x, and the indexes a ^ /3 7 = 1, 2, 3. The functions 

fa{y) satisfy the following boundary conditions fa{y) 0, when y —> 00 . For a system including particles with spin 
and isospin, the number of equations in the set depends on the number of terms of the Faddeev components ()20l) in 
the spin-isospin basis. For the K~pp system considered within the framework of the s-wave approach, we have the 
set of three equations corresponding to (fTO . The set of eigenvalues of pair subsystems e^, i = 1,2,3 includes €2 as 
the bound state energy of the K~p pair in singlet isospin state assumed to be the A(1405) hyperon. 


III. RESULTS AND DISCUSSION 

Results of our calculations for the K~pp cluster are presented in this section. For the calculations of the binding 
energy and the width with the methods of HH and DFE we use for the NN interaction MT, T and AV14 potentials, 
while for the KN interaction we use the single-channel phenomenological AY KN potential with the range parameter 
b = 0.66 fm and the energy dependent effective HW potential with the range parameter 6=0.47 fm. Such an approach 
allowed us to examine how the K~pp cluster’s structure depends on different choices of the KN interactions for the 
same NN potential, as well as dependence from different choices of the NN interaction for the same KN interaction. 
This enable us to understand the sensitivity of the system to the input interactions. By solving the system of equations 
(El one finds the binding energy of the K pp cluster. The results of the calculations with the method of HH strongly 
depend on the number of terms in the expansion of the wave function (nnj. The convergences of binding energy 
calculations for the ground state of the K~pp system as a function of the grand angular momentum y, are shown in 
Fig. [2] for different input potentials. Reasonable convergence is reached for pmax = 10 and we limit our considerations 
to this value. As is seen from Fig. [5] much faster convergence occurs for considered NN potentials when the AY 
potential is chosen for the KN interaction, while for the HW KN interaction the binding energy calculations converge 
more slowly. The accuracy reached is about 0.2 MeV and our expectation is that consideration of the higher values 
of the grand angular momentum y would not dramatically change the binding energy of the K~pp cluster. 

The results for the binding energy and the width of the K~pp system using the method of HH for different KN 
and NN interactions are presented in Table HI Ej^-p is the two-body energy when only the K~p pair are interacting, 
while the interaction between the two protons is neglected. The HH method allows one to obtain the wave function of 
the K~pp cluster. By solving the system of equations El one finds the binding energy as well as the corresponding 
hyperradial functions. The latter allows one to construct the wave function d)* for the K~pp system. Using the wave 
function, the width of the bound state can be evaluated in a perturbative way from the imaginary part of the KN 
interaction as F = —2 (41 \lTtr{yK-p(r 12 ) + Uif-p(r’i 3 )) | T). Since |lmyR--p(r-)| |ReVR--p(r)|, this is a reasonable 
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FIG. 2: The convergence of the ground state energies E of the K pp system as a function of pmax for different NN potentials 
and K~ N interactions. 

TABLE I: The binding energy B and width T for the K~pp system calculated in the framework of the method of HH in the 
momentum representation for different interactions. NN potentials: MT and T [^. KN interactions: AY @] and HW 
[^ . Ef^-p is two-body energy of the K~p. 



MT+AY T-f AY 

MT-bHW 

T-bHW 

B, MeV 

46.5 46.3 

20.5 

20.6 

r, MeV 

84.3 74.5 

48.1 

49.5 

MeV 

29.9 

10.9 


approximation for the width. As it is stated in Gal’s review , as well as demonstrated in the recent calculations of 
the width for the K~pp system using a coupled-channel complex scaling method with Feshbach projection, this 
is a good approximation. For an approximate evaluation of the width the imaginary part of the complex potential 
has often been treated perturbatively by many authors [III, [H, Ibll . . Recently this approach was used in 

Ref. [13 for few-body calculations of ? 7 -nuclear quasibound states. We use the latter expression to find the width and 
the calculated values of the width for different KN and NN interactions are listed in Table HI The analysis of the 
results presented in Table U shows that, for the chiral HW KN interaction, the binding energy is less than half that 
of the phenomenological AY interaction and the HW interaction also leads to a much smaller width. However, the 
binding energy and the width are not sensitive to the form of NN potentials for the same type oi KN interaction. 
The binding energy turns out to be about 1 MeV, while the width ranges from ^75 to ^84 MeV for the AY potential 
and turns out to be within about 1.5 MeV for the HW KN interaction. 

We will now discuss the results when we employ the Faddeev technique. The results for the binding energy of 
the K~pp system for the orbital momentum configurations of {pp)K~ and {K~p)p rearrangements are presented in 
Table m The calculations are performed using the AV14 NN potential for the AY and HW KN interactions. As 
can be seen from Table El the orbital contributions with {l,X) are equal (2,2) and (4,4) are small enough and the 
s—wave consideration is very reasonable. The numerical convergence for the binding energy calculations is fast with 
increasing model space. The visible contribution that provides the bound state of the system comes from the p-wave 
of the pp pair. For the AV14 potential, we also calculated the binding energy using the averaged potential approach 
as well as the cluster approach. The results with the averaged potential (IT^ are 33.7 MeV and 12.5 MeV for the 
AY and HW KN interactions, respectively. We used dominant clustering of the K~pp system in the configuration 
A(1405) -bp to calculate the binding energy within the cluster approach for the DFE. The Faddeev components are 
decomposed as a product of the eigenfunctions of the Hamiltonians for the two particle subsystems and the function 
of the relative motion of the pairs and the third particle. Our results for the functions associated with the relative 
motion are shown in Fig. [3] for the nucleon-nucleon MT potential and the AY and HW effective KN interactions. 
One can evaluate the dominant contribution of the (iT“p)®™®*®‘ +p configuration to the total wave function of the 
system, which is represented as a sum of the Faddeev components. The calculations for other NN potentials have 
shown the same behavior and relative contributions of the functions associated with relative motion. In the cluster 
approach, the binding energy increases for both K~p interactions and are 49.5 MeV and 17.7 MeV, respectively. 

Comparing the binding energies obtained by using the averaged potential and the cluster approaches with the s- 
wave DFE calculations, one can conclude that the averaged potential approach gives the better result for the relatively 
weaker HW potential. At the same time, the cluster approach is more accurate for the stronger AY potential. We also 
calculated the proton separation energy within the averaged potential, the cluster and the s-wave DFE approaches. 
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TABLE II: Orbital momentum configurations of {pp)K~ and {K~p)p rearrangements and the ground state energy of the K~pp 
system calculated with the HW and AY KN potentials. The NN interaction is given by the AV14 potential [s^. I is the 
orbital momentum of the pair of particles, and A is the orbital momentum associated with the relative motion of the third 
particle with respect of the center of mass of the pair. 

(pp)K- _ {K-p)p HW AY 

[I, A) (0,0) (0,0) -21.15 -46.97 

{I, A) (0,0) (2,2) (0,0) -21.54 -47.33 

jl, A) (0,0) (2,2) (4,4) (0,0) -21.56 -47.34 


0.15 

0.12 
_ 0.09 

3 

nj 

X 0.06 
0.03 

0.00 

-0.03 

0 5 10 15 0 5 10 15 

y (ftn) y (fm) 

FIG. 3: The function of relative motion within the cluster approach for the K~pp system: the dot-dashed curve shows 
calculation for the configuration K~ + {pp)\ the solid curve shows calculation for the configuration {K~p) +p with the singlet 
isospin state of the pair {K~p); the dashed curve shows calculation for the configuration {K~p) +p with the triplet isospin 
state of the pair {K~p). The configurations are related to the set of Eqs. (1211) . The calculations are performed with the 
nucleon-nucleon MT interaction for a) the AY effective KN potential and b) the HW effective KN potential. 



The results are presented in Table uni 

The summary of the results for the binding energy of the K~pp system in the framework of the method of DFE 
for all considered NN and KN interactions are presented in Table El The analysis of the calculations presented in 
Table HVl shows that the AY potential as the KN interaction input fall into the 46 - 47 MeV range for the binding 
energy of the K~pp cluster, while the chiral HW KN potential gives about 20.4 - 21.6 MeV for the binding energy. 
Now one can address the theoretical discrepancies in the binding energy for the K~pp system related to the different 
NN and KN interactions. The comparison of the results of calculations presented in Tables HI and IIVI for the binding 
energy for the K~pp system obtained by both methods are in reasonable agreement. The ground state energy is 
not sensitive to the NN interaction. However, it shows very strong dependence on the kaon-nucleon potential. The 
energy of the ground state, as well as the width calculated for the e nerg y-independent KN interaction Q are more 
than twice bigger than for the energy-dependent chiral KN potential [^. Therefore, the highest binding energies are 
those that are obtained based on the phenomenological AY potential. Discrepancies obtained for the binding energy 


TABLE III: The proton separation energy calculated for the MT NN potential and the HW and AY KN interactions. 

Energies are given in nnit of MeV. 


KN 

Averaged potential 

Cluster 

DFE 

interaction 

approach 

approach 


AY 

24.1 

19.2 

15.7 

HW 

11.2 

6.5 

9.3 
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TABLE IV: The binding energy B for the K pp system calculated in the framework of the DFE method for different NN and 
KN interactions. NN potentials: AV14 [s^, MT and T [hJl. KN interactions: AY and HW [s^. a is the two-body 
scattering length and is the root-mean-squared distance in the K~p system. is the two-body energy of K~p. 


KN 

a, fm 

, fm 

MeV 

NN 

B, MeV 

AY 

1.88 

1.35 

30.26 

AV14 

47.3 





MT 

46.0 





T 

46.3 

HW 

2.68 

1.94 

11.16 

AV14 

21.6 





MT 

20.4 





T 

20.6 


TABLE V: Summary of the theoretical studies for the Kpp cluster. 


Method 

B{K-pp) Width, r 

KN 

References 


MeV 

MeV 



Variational 

48 

61 

AY 

[8], [9], [1^ 

Methods 

20±3 

40-70 

Chiral model 

[33], [3^ 


40-80 

40-85 

Separable 

m 


20-35 

20-65 

Chiral model 

[37] 


124 

12 

AY 

m 

Methods of 

47-70 

90-100 

Separable En. Indep. 

i 

i 

i 

Faddeev 

~32 

-50-65 

Separable En. Dep. 

m 

equations 

45-95 

45-80 

Separable En. Indep. 

[42], [4^ 


9-16 

34-40 

Separable En. Dep. 

m 


30-40 

50-80 


m - [4^ 


-52 


Separable En. Indep. 

[49] 

Methods of HH 

-16 

-41 

Chiral model 

m 


15-17 

36-43 

Chiral model 



40-48 

75-96 

AY 



using the same potentials but different methods, - the method of HH in momentum representation and the method 
of Faddeev equations in configuration space,- are mostly related to a problem of an equivalent representation of the 
potentials in momentum and configuration spaces. 

Now let us compare our results with those obtained using variational methods, the method of Faddeev equations in 
the momentum representation and the method of hyperspherical functions in configuration space. The summary of the 
results for the binding energy and the width with different theoretical approaches and models, obtained for different 
kinds oi NN and KN interactions, are presented in Table E We start the comparison with variational calculations. 
In Ref. the variational calculations were carried out by constructing a wave function for the three-body system 
with correlation functions for each of the constituent pairs on the basis of multiple scattering theory. The binding 
energy and the width for the K~pp cluster was calculated by employing the AY potential as the KN interaction 
and the bare Tamagaki GSRS potential as the NN interaction. Thus, the authors of [T^ used the same s-wave 
interactions, and this allows us to compare our results with the previous one. The binding energy found with the DFE 
and HH methods are in good agreement with the one obtained with the variational method. This is a good sign that 
the binding energy does not depend significantly on the method of calculation. The prediction in Ref. [s^ for the 
binding energies of Itot = 1/2, Inn = 1 states given by the s-wave interactions and described by multiple scattering 
in the single KN channel falls into the 40 - 80 MeV range with the parametrization for the Arg onne AV18 potential 
from Ref. . These results are consistent with our calculations using the AV14 potential [HJ and the results from 
[a. The differences within this range are mostly due to a different KN input and possibly slightly due to the NN 
input. 

Different variational approaches used in Refs. [1^ and [s^l are of comparable quality in their high degree of 
consistency. Dote, Hyodo and Weise [s^l employed several versions of energy-dependent effective KN interactions 
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derived from chiral SU(3) dynamics together with the realistic Argonne 18 NN potential They found that the 
antikaonic dibaryon K~pp is not deeply bound and obtained 20±3 Mev for the binding energy. Our calculations in 
the framework of the DFE and HH methods when we employ the effective energy-dependent chiral-theory based HW 
potential for KN interaction and different NN interactions, as inputs, also predict a shallowly bound K~pp cluster. 
This is consistent with results from [s^, [sj and recent calculations (s^ . 

The first calculations with Faddeev equations in Alt-Grassberger-Sandhas form for the three-body system with 
coupled KNN and ttYiN channels were performed in Re fs. 13^ . [40l | for the study of the KNN system with separable 
two-body potentials yield large bindings. In Refs. l40j| authors obtained 55-70 MeV, and 95-110 MeV for the 
binding energy and the width, respectively. With a similar approach, Ikeda and Sato [4l| calculated B ~80 Mev and 
r ~73 MeV. Later, two of the authors of j4l| repeated their calculation in [H, using two models with the energy- 
independent and energy-dependent potentials for the s—wave KN interaction, and their calculations yield smaller 
values for the binding energy 44-58 MeV and width 34-40 MeV [d^. Recently in Ref. [s^, the Faddeev calculations 
for the KNN quasi-bound state with the two phenomenological and the energy-dependent chirally motivated models 
of the KN interaction were carried out. The binding energy for the K~pp cluster obtained was 32 MeV with the 
chirally motivated models and 47 - 54 MeV with the phenomenological KN potentials. The width is about 50 - 
65 MeV. Therefore, we can conclude that the Faddeev calculations for the energy-independent models for the KN 


interaction 

interaction 



1^ . Id^ predict a deeper binding energy than that of the energy-dependent description of the KN 
. Our calculations obtained with both methods confirm that the effective KN interaction derived 


from chiral SU(3) dynamics yields a shallowly bound K pp cluster, while the phenomenological energy independent 


AY potential predicts much deeper binding energy 


our results with the Faddeev calculations _ 
to describe the KN and NN interactions. 

The fixed center approximation (FCA) to the Faddeev equations [4 


for all considered NN interactions. A more detailed comparison of 
y, l50j is not easy because these calculations use separable potentials 


48j| leads to the binding energy of the K 


system with S' = 0 30-40 MeV and the width 50-80 MeV that includes the recent improvement of the method [47|, 14^ 
by including the charge exchange mechanisms in the K~ rescattering and absorption which have been ignored in 
previous works within this approximation. Due to the absorption of K~ by two nucleons, the width of the bound 
K~pp cluster is increased by about 30 MeV [i^. Although the FCA to the Faddeev equation makes a static picture of 
the two nucleon and does not consider the recoil of the spectator nucleon, this approach provides reasonable values for 
the binding energy and width, as seen from Table IVl and qualitatively corroborates findings done with other methods 
which are technically much more involved. 

The study of K~pp with 5 = 0 was performed using the method of hyperspherical functions in configuration space 
51 1 and in the momentum representation . In Ref. (Fij i , the binding energies and widths of the three-body K~pp 
cluster are calculated using the realistic AV4 NN potential and a subthreshold energy-dependent chiral KN 
interaction derived with a chiral model [30l| . The results are in good agreement with previous K~pp calculations 
with an energy-dependent chiral KN interaction as input. Calculations of the binding energy and the width 
with the method of HH in the momentum representation |52| reproduced the results of Ref. (5ll |. and we assume 
the slight difference arised from the difference between using AV4 and AV18 potentials and perhaps is related to the 
conversion of the hyperspherical expansion. However, most importantly, our results and those of the previous HH 
calculations (Fll . [F^ j are in good agreement and support the conclusion that the key role in binding the K~pp system 
is played by the KN interaction and the KN potential obtained based on chiral SU(3) dynamics leads to binding 
energies of relatively low values. Also, these calculations show that the binding energies are small and the widths 
are more than twice that of the binding energies. Therefore, all model calculations with the KN interaction derived 
based on the chiral unitary approach predict a shallow binding state with very large width for the K~pp cluster. This 
is a precautionary indication that it may be difficult to experimentally observe the K~pp cluster. 


Conclusion 


Within the framework of a potential model for the kaonic cluster K~pp, we perform nonrelativistic three-body 
calculations using two methods: the method of hyperspherical harmonics in the momentum representation and the 
method of Faddeev equations in configuration space. We examine how the quasi-bound state of the K~pp cluster 
depends on different choices of the KN and NN interactions. Our consideration includes the realistic AV14 [F^, the 
semi-realistic Malfliet and Tjon MT-I-HI [F^ and Tamagaki C3RS m potentials as input for the NN interaction 
and we employ the phenomenological AY potential [ 1 , [3 HW potential [F^ based on chiral SU(3) dynamics 
as input for the KN interaction. The results for the binding energy of the K~pp system obtained by the method 
of HH and the DFE method are in reasonable agreement. A discrepancy obtained is related to a problem of the 
equivalent representation of the potentials in momentum and configuration spaces. For all types of considered NN 
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interactions, both methods predict deeply bound states for the AY KN interaction and a relatively shallowly bound 
K~pp cluster for the effective KN interactions derived from chiral SU(3) dynamics. Moreover, the K~pp cluster is the 
most strongly quasi-bound three-body system. The results of our calculations within the DFE and HH methods show 
that the binding energy of the K~pp system depends entirely on the ansatz for the KN interaction and substantially 
changes when we use a phenomenologically constructed KN potential [l^ and a KN potential obtained within the 
framework of the chiral unitary approach pill. Perhaps the ambiguity of the KN interaction stems from an 
accuracy of description of experimental data for the energy shift in the kaonic hydrogen atom and K~p scattering 
by the KN potentials. Related to the sensitivity of the binding energy to the details of the NN potentials Ref. 
[301 mentioned that as long as the K~pp system is only weakly bound, the dependence on different types of NN 
interactions is weak. In fact, our study confirms this conclusion and, moreover, shows that the dependence on different 
types of NN interactions is also weak if the K~pp system is strongly bound. Using the formalism of the Faddeev 
equations, we may separate channel configurations of the total wave function of the system. We have shown that the 
configuration {K~p) + p with the singlet isospin state of the pair K~p dominates in the system due to strong K~p 
interaction. We have seen that, on the background of the KN potential, different NN interactions weakly effluent to 
the bound state energy. It has to be noted that the s-wave models, with the simplifications taken, provide reasonable 
description of the system. The analysis of data presented in Table |V| shows that the width is always larger than the 
binding energy. Particularly for some calculations, the width is more than twice as much as the binding energy. 
Thus, we are facing a situation in which the states have a much larger width than binding energy, which makes the 
experimental observation challenging [4^ and it may be hard to identify the resonance. However, the continuation of 
the experimental search for the quasi-bound kaonic cluster still remains important. 
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